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Abstract We calculate in detail the expected properties of low redshift DLAs under the 
assumption that they arise in the gaseous disks of galaxies like those in the 2: ~ 
population. A sample of 355 nearby galaxies is analysed, for which high quality 
Hi 21 -cm emission line maps are available as part of an extensive survey with 
the Westerbork telescope (WHISP). We find that expected luminosities, impact 
parameters between quasars and DLA host galaxies, and metal abundances are 
in good agreement with the observed properties of DLAs and DLA galaxies. 
The measured redshift number density of 2: = gas above the DLA limit is 
dN/dz = 0.045 ± 0.006, which compared to higher z measurements implies 
that there is no evolution in the comoving density of DLAs along a line of sight 
between 2 ~ 1.5 and z = 0, and a decrease of only a factor of two from 2 ~ 4 
to the present time. We conclude that the local galaxy population can explain all 
properties of low redshift DLAs. 
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1. Introduction 

The range of Hi column densities typically seen in routine 21 -cm emis- 
sion line observations of the neutral gas disks in nearby galaxies is very sim- 
ilar to those that characterise the Damped Lyman-a Systems or DLAs with 
A^Hi > 2 X 10^'^ cm^^. An attractive experiment would therefore be to map 
the Hi gas of DLA absorbing systems in 21 -cm emission, and measure the 
DLAs' total gas mass, the extent of the gas disks and their dynamics. This 
would provide a direct observational link between DLAs and local galaxies, 
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but unfortunately such studies are impossible with present technology (see 
e.g., Kanekar et al. 2001). The transition probability of the hyperfine split- 
ting that causes the 21-cm Hne is extremely small, resulting in a weak line that 
can only be observed in emission in the very local {z < 0.2) universe, with 
present technology. On the other hand, the identification of DLAs as absorbers 
in background QSO spectra is, to first order, not distance dependent because 
the detection efficiency depends mostly on the brightness of the background 
source, not on the redshift of the absorber itself. In fact, the lowest redshift 
(z < 1.7) Lyman-a absorbers cannot be observed from the ground because the 
Earth's atmosphere is opaque to the UV wavelength range in which these are to 
be found. Furthermore, due to the expansion of the universe the redshift num- 
ber density of DLAs decreases rapidly toward lower redshifts. Consequently, 
there are not many DLAs known whose 21-cm emission would be within the 
reach of present-day radio telescopes. 

So, we are left with a wealth of information on the cold gas properties in 
local galaxies, which has been collected over the last half century, and several 
hundreds DLA absorption profiles at intermediate and high redshift, but little 
possibility to bridge these two sets of information. Obviously, most observers 
resort to the optical wavelengths to study DLAs but attempts to directly image 
their host galaxies have been notably unsuccessful (see e.g.. Warren et al. 2001 
and M0ller et al. 2002 for reviews). A few positive identifications do exist, 
mostly the result of HST imaging. 

Although the absolute number of DLAs at low z is small, the success rate for 
finding low-z host galaxies is better for obvious reasons: the host galaxies are 
expected to be brighter and the separation on the sky between the bright QSO 
and the DLA galaxy is hkely larger Early surveys for low-z DLA host galax- 
ies consisted of broad band imaging and lacked spectroscopic follow-up (e.g., 
Le Brun et al.l997). Later studies aimed at measuring redshifts to determine 
the association of optically identified galaxies with DLAs, either spectroscop- 
ically (e.g., Rao et al. 2003), or using photometric redshifts (Chen & Lanzetta 
2003). All together, there are now ~ 20 DLA galaxies known at z < 1. The 
galaxies span a wide range in galaxy properties, ranging from inconspicuous 
LSB dwarfs to giant spirals and even early type galaxies. Obviously, it is not 
just the luminous, high surface brightness spiral galaxies that contribute to the 
H I cross section above the DLA threshold. As explained above, we cannot 
study these galaxies in the 21-cm line on a case-by-case basis, but we can do a 
study of a statistical nature to see if the properties of DLAs and DLA galaxies 
agree with our knowledge of H i in the local universe. 
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2. 



140,000 "DLAs" at 2 = 



Blind 21 -cm emission line surveys in the local universe with single dish 
radio telescopes such as Parkes or Arecibo have resulted in an accurate mea- 
surement of Q.Bi{z = 0), which can be used as a reference point for higher 
redshift DLA studies, fini is simply calculated by integrating over the Hi 
mass function of galaxies, which is measured with surveys such as HIPASS 
(Zwaan et al. 2005a). However, due to the large beam widths of the singe dish 
instruments, these surveys at best only barely resolve the detected galaxies and 
are therefore not very useful in constraining the column density distribution 
function of z H I. Hence, for this purpose we use the high resolution 21- 
cm maps of a large sample of local galaxies that have been observed with the 
Westerbork Synthesis Radio Telescope. This sample is known as WHISP (van 
der Hulst et al. 2001) and consists of 355 galaxies spanning a large range in H i 
mass and optical luminosity. The total number of independent column density 
measurements above the DLA limit is ^ 140, 000, which implies that the data 
volume of our present study is the equivalent of ~ 140, 000 DLAs at z = 0! 

Each galaxy in the sample is weighted according to the H i mass function 
of galaxies. We can now calculate the column density distribution function 



where S(A'hi, Afni) is the area function that describes for galaxies with Hi 
mass Mhi the area in Mpc~^ corresponding to a column density in the range 
A^Hi to A^Hi + dNm, and $(Mhi) is the H i mass function. c/Hq converts the 
number of systems per Mpc to that per unit redshift. 

Figure 1 shows the resulting /(-/Vhi) on the left, and the derived Hi mass 
density per decade of N-ai on the right. For comparison with higher redshift 
observations, we also plot the results from two other studies. The Peroux 
(2005) measurements of /(A^hi) below the DLA limit are the result of their 
new UVES survey for "sub-DLAs". The intermediate redshift points from 
Rao et al. (2005) are based on Mg ll-selected DLA systems. The surprising re- 
sult from this figure is that there appears to be only very mild evolution in the 
intersection cross section of H i from redshift z ~ 5 to the present. From this 
figure we can determine the redshift number density of log iVni > 20.3 gas 
and find that dN/dz = 0.045 ± 0.006, in good agreement with earlier mea- 
surements at z = 0. Compared to the most recent measurements of dN/ dz at 
intermediate and high z, this implies that the comoving number density (or the 
"space density times cross section") of DLAs does not evolve after z ~ 1.5. In 
other words, the local galaxy population explains the incidence rate of low and 
intermediate z DLAs and there is no need for a population of hidden very low 
surface brightness (LSB) galaxies or isolated H i clouds (dark galaxies). 
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Figure 1. Left: The Hi column density distribution function f{Nm)- The solid dots are 
from our analysis of 21-cm maps of local galaxies, the other points are from various surveys at 
higher redshifts as indicated in the legend. All calculations are based on a Qm = 0.3, Qa = 0.7 
cosmology. Right: The contribution of different column densities (per decade) to the integral 
H I mass density. Note the very weak evolution in flm from high z to the present. 



The right hand panel shows that at z = most of the H i atoms are in col- 
umn densities around A'hi = 10^^ cm^^. This also seems to be the case at 
higher redshifts, although the distribution might flatten somewhat. The one 
point that clearly deviates is the highest A'^hi point from Rao et al. (2005) at 
log A'^Hi = 21.65. The figure very clearly demonstrates that this point domi- 
nates the 0,111 measurement at intermediate redshifts. It is therefore important 
to understand whether the Mg ii-based results really indicate that high column 
densities (log A'^hi 21.65) are rare at high and low redshift, but much more 
ubiquitous at intermediate redshifts, or whether the Mg ii selection introduces 
currently unidentified biases. 

3. Expected properties of low-z DLAs 

Now that we have accurate cross section measurement of all galaxies in our 
sample, and know what the space density of our galaxies is, we can calculate 
the cross-section weighted probability distribution functions of various galaxy 
parameters. Figure 2 shows two examples. The left panel shows the i?-band 
absolute magnitude distribution of cross-section selected galaxies above four 
different H i column density cut-offs. 87% of the DLA cross-section appears to 
be in galaxies that are fainter than L^, and 45% is in galaxies with L < L^,/10. 
These numbers agree very well with the luminosity distribution of z < 1 DLA 
host galaxies. Taking into account the non-detections of DLA host galaxies and 
assuming that these are <C L*, we find that 80% of the z < 1 DLA galaxies are 
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Figure 2. Left: The expected distribution of B-band absolute magnitudes of z = high H I 
column density systems. The lines plus errorbars show the product of cross sectional area and 
space density, which translates to the number of expected absorbers per Mpc per magnitude. 
The right axis shows the corresponding number of absorbers per unit redshift dN/dz. The 
different lines correspond to different column density limits, as indicated by the labels. The 
thick line corresponds to the classical DLA limit of log A'^ei > 20.3. Right: The probability 
distribution of oxygen abundance [O/H] of H I absorbers. The solid lines refer to the approach 
of assuming fixed [O/H] gradients and the the dashed lines refer to varying gradients (see text). 



sub-L*. The median absolute magnitude of a z = DLA galaxy is expected 
to be Mb = —18.1 (~ -L*/7). The conclusion to draw from this is that we 
should not be surprised to find that identifying DLA host galaxies is difficult. 
Most of them (some 87%) are expected to be sub-L* and many are dwarfs. 

Using similar techniques, we find that the expected median impact param- 
eter of log A^Hi > 20.3 systems is 7.8 kpc, whereas the median impact pa- 
rameter of identified z < \ DLA galaxies is 8.3 kpc. Assuming no evolution 
in the properties of galaxies' gas disks, these numbers imply that 37% of the 
impact parameters are expected to be less than 1" for systems at z = 0.5. This 
illustrates that very high spatial resolution imaging programs are required to 
successfully identify a typical DLA galaxy at intermediate redshifts. 

The right panel in Figure 2 shows the probability distribution of oxygen 
abundance in z = DLAs. We constructed this diagram by assigning to every 
Hi pixel in our 21 -cm maps an oxygen abundance, based on the assumption 
that the galaxies in our sample follow the local metallicity-luminosity {Z — L) 
relation (e.g. Gamett 2002), and that each disk shows an abundance gradi- 
ent of [O/H] of — 0.09 dexkpc"^ (e-g-, Ferguson et al. 1998) along the major 
axis. The solid lines correspond to these assumptions, the dotted lines are 
for varying metallicity gradients in disks of different absolute brightness. The 
main conclusion is that the metallicity distribution for H i column densities 
log A^Hi > 20.3 peaks around [0/H]=— 1 to —0.7, much lower than the mean 



6 



value of an galaxy of [0/H]« 0. The reason for this being that 1) much 
of the DLA cross section is in sub-L* galaxies, which mostly have sub-solar 
metallicities, and 2) for the more luminous, larger galaxies, the highest inter- 
ception probability is at larger impact parameters from the centre, where the 
metallicity is lower. Interestingly, this number is very close to the z = extrap- 
olations of metallicity measurements in DLAs at higher z from Prochaska et al. 
(2003) and Kulkarni et al. (2005). For the mean mass-weighted metallicity of 
H I gas with log Nm > 20.3 at z = we find the value of log Z = -0.35±0.2, 
also consistent with the z = extrapolation of iVni-weighted metallicities in 
DLAs, although we note this extrapolation has large uncertainties given the 
poor statistics from DLAs at z < 1.5. These results are in good agreement 
with the hypothesis that DLAs arise in the H I disks of galaxies. 

4. Conclusions 

The local galaxy population can explain the incidence rate and metallicities 
of DLAs, the luminosities of their host galaxies, and the impact parameters 
between centres of host galaxies and the background QSOs. 

This work is presented in much more detail in a forthcoming paper (Zwaan 
et al. 2005b). 
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